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ABSTRACT 

We report on the discovery of a new quantitative relationship between X-ray gas morphology and 
radio and X-ray AGN luminosities in normal elliptical galaxies. This is the second paper in a series 
using data on 54 objects from the Chandra public archive and builds on the findings of Paper I, which 
demonstrated that hydrostatic equilibrium in elliptical galaxies holds, at best, only approximately 
and that the shape of the X-ray isophotcs is unrelated to the shape of the gravitational potential. 
Instead, the gas is almost always asymmetrically disturbed. In this paper, we quantify the amount of 
asymmetry and study its correlation with other galaxy properties. We also determine radio powers 
and derive X-ray AGN luminosities for our galaxy sample. We find that the amount of asymmetry 
in the gas is correlated with both measures of AGN activity, in the sense that the hot gas is more 
disturbed in galaxies with higher radio and X-ray AGN luminosities. We find no evidence that galaxy 
density has significant effects on gas morphology. We do however find evidence for a correlation 
between gas asymmetry and the presence of hot ambient gas, which we interpret as a signature of 
hydrodynamic interactions with an external ambient medium. Surprisingly, the AGN-morphology 
connection persists all the way down to the weakest AGN luminosities in rather X-ray faint galaxies. 
This is strong morphological evidence that supports the general importance of AGN feedback, even 
in normal elliptical galaxies. 

Subject headings: galaxies: cooling flows — galaxies: elliptical and lenticular, cD — galaxies: ISM — X- 
rays: galaxies — X-rays: ISM 
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1. INTRODUCTION 

In the first paper of this series (Dichl & Statler 2007, 
hereafter Paper I) we conducted a morphological analysis 
of the hot interstellar medium (ISM) of normal elliptical 
galaxies, derived from archived Chandra observations of 
54 objects. This paper introduced a new technique to 
isolate the diffuse hot gas emission from the contami- 
nating effects of unresolved point sources, and presented 
a gallery of adaptivcly binned gas-only images. By fit- 
ting elliptical isophotes to these images, we demonstrated 
that the apparent flattening of the X-ray gas is com- 
pletely uncorrelated with, and often significantly larger 
than, that of the starlight. The absence of any correla- 
tion, at radii inside the optical effective radius, implies 
that the hot gas cannot be in perfect hydrostatic equilib- 
rium, unless the potential is dominated by a dark mat- 
ter component completely unrelated to the stellar mass 
distribution. This is unlikely, since within the effective 
radius gravitational potentials are generally expected to 
be stellar-mass dominated (e.g. Mamon & Lokas 2005; 
Humphrey et al. 2006). Consequently, efforts to measure 
the shapes of dark matter halos using X-ray gas isophotes 
may be fruitless in most cases (e.g. Buote & Canizares 
1994; Buote et al. 2002), and radial mass profiles derived 
from X-ray data (Humphrey et al. 2006; Fukazawa et al. 
2006) could be in error by factors of order unity. 

This paper takes up the question left unaddressed by 
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Paper I: what is the origin of the gas fiattcning, if not 
gravitational? In Paper I, we already excluded rotational 
fiattening as the dominant cause. A qualitative analysis 
of gas morphologies showed that they are often asym- 
metric and disturbed. In this paper we will quantify 
this asymmetry, argue that it and the isophotal cUip- 
ticity refiect the same underlying disturbance, and show 
strong evidence that the root causes are the central ac- 
tive nucleus — even in those systems with very weak ac- 
tive galactic nuclei (AGN) — and interactions with ambi- 
ent external gas. 

Asymmetries are well known in X-ray studies of galax- 
ies, groups, and clusters. Chandra and XMM-Newton ob- 
servations of clusters of galaxies reveal strong deviations 
from smooth surface brightness distributions. Many of 
these features are believed to be due to cold fronts, 
caused by infalling substructure into the cluster center 
(e.g. Chatzikos et al. 2006), or gas sloshing in the grav- 
itational potential (e.g. Ascasibar & Markevitch 2006). 
In other cases, depressions in the surface brightness dis- 
tribution are found to be coincident with extended radio 
emission (e.g. Nulsen et al. 2005). These depressions, or 
"bubbles" , are generally believed to be pockets of low- 
density extremely hot plasma, inflated by jets powered 
by the central AGN (see McNamara & Nulsen 2007, for a 
recent review). However, there are many cases where the 
depressions have no obvious radio counterparts, which 
has brought them the nickname "ghost cavities" (McNa- 
mara et al. 2001). These cavities are now interpreted to 
be relic bubbles that have detached from the radio source 
and are buoyantly rising radially outward in the cluster 
gas. 

Unexpectedly, Chandra and XMM observations of nor- 
mal elliptical galaxies have also revealed a wealth of 
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highly disturbed gas morphologies with qualitatively 
similar properties to clusters. Many individual obser- 
vations have been analyzed in detail and a variety of 
different explanations have been proposed. Sharp, one- 
sided drops in surface brightness are usually interpreted 
as signs of ram-pressure, distorting the gas in ellipti- 
cal galaxies as they move through the ambient intr- 
acluster or intragroup medium, as in NGC 4472 (Ir- 
win & Sarazin 1996; BiUer et al. 2004) and NGC 1404 
(Machacek et al. 2004). To explain the very asymmet- 
ric emission in NGC 7618, Kraft et al. (2006) even ar- 
gue for a major group-group merger. Other observations 
are not interpreted as results of environmental effects, 
but rather as signatures of the central AGN. For a few 
cases, such as NGC 4374 (Finoguenov & Jones 2001) or 
NGC 4472 (Biller et al. 2004), the association with the 
AGN is indicated clearly by the correspondence between 
the radio source and depressions in the X-ray gas dis- 
tribution. While these cavities are generally surrounded 
by cool rims (McNamara & Nulsen 2007, and references 
therein), two low-power radio galaxies (Centaurus A and 
NGC 3801) have recently been identified with shocks sur- 
rounding the radio lobes (Kraft et al. 2003, 2007; Croston 
et al. 2007). In other objects, one can directly detect the 
X-ray counterpart of the radio jet (e.g. Sambruna et al. 
2004; Harris et al. 2002; Hardcastle et al. 2003). However, 
most cases are less clear-cut. For example, the origin of 
the shock-like features seen in NGC 1553 (Blanton et al. 
2001) or NGC 4636 (Jones et al. 2002) is still a mystery 
The morphologies in other galaxies qualitatively resem- 
ble ghost cavities very close to the core (e.g. NGC 5044, 
Buote et al. 2003), without evidence for an active radio 
source. 

These individual cases each come with their own analy- 
sis techniques and interpretations, making general state- 
ments about elliptical galaxies problematic. Detailed 
analyses are generally limited to the X-ray brightest 
galaxies with clearly identifiable features, which are often 
dismissed as special cases. However, our morphological 
survey of normal elliptical galaxies in Paper I shows that 
disturbances and asymmetries in the gas are actually the 
norm rather than the exception, even for relatively X-ray 
faint galaxies. Most galaxies in the Chandra archive lack 
the signal, and therefore the contrast, to reliably identify 
cavity or shock features individually. In this paper, we 
introduce a statistical measure of asymmetry (the asym- 
metry index, rj) that allows us to treat all objects on an 
equal footing. We will show that the asymmetry index is 
strongly correlated with two different measures of nuclear 
activity, and with the presence of ambient hot gas. 

The analysis in this paper begins with the elliptical 
isophotal fits of Paper I, and is organized as follows. In 
§2 we review the techniques of Paper I and describe the 
sample. We then calculate smooth functional fits to the 
isophotal profiles, use these fits to generate symmetric 
surface brightness models, and use the residuals from 
these models to define the asymmetry index rj for each 
galaxy. We show by simulations that 77 is a monotonic 
measure of the sort of asymmetry expected from multiple 
cavities, and that it is independent of non-morphological 
parameters. We describe our techniques to extract the 
central AGN X-ray luminosity and to determine radio 
power and environment. In §3 we show that the asym- 
metry and isophotal ellipticity are correlated, and thus 



probably measure the same underlying disturbance. We 
show that neither of these quantities depends on neighbor 
galaxy density, but that they are correlated with mea- 
sures of AGN power and the presence of ambient gas. 
The consequences for AGN feedback in galaxies and in- 
teractions with the ambient medium are discussed in §4, 
before we conclude with a reiteration of our major results 
in §5. 

2. DATA ANALYSIS 
2.1. Summary of Paper I and Preliminary Analysis 

Our full sample, as described in Paper I, consists of 
54 early-type galaxies observed with the the ACIS-S 
instrument on the Chandra satellite during cycles 1- 
4. The data have been homogeneously reprocessed to 
avoid problems due to changes in the standard calibra- 
tion pipeline over time. We will briefly summarize the 
generation of those data products essential for the un- 
derstanding of this paper and refer the interested reader 
back to Paper I for more details. 

In Paper I, we presented a new method to isolate the 
hot gas emission in elliptical galaxies from the contami- 
nation of unresolved point sources. This method is based 
on the fact that point sources and hot gas contribute dif- 
ferently to soft (0.3 - 1.2 keV) and hard (1.2 - 5keV) 
bands, due to their intrinsically different spectra. Thus, 
by subtracting a properly scaled version of the hard band 
image from the soft band after removing resolved sources, 
we are able to remove the contribution of the harder 
unresolved point sources and isolate the soft gas emis- 
sion self-consistcntly, correcting for the amount of gas 
flux subtracted from the hard band as well. All im- 
ages are photon-flux calibrated, corrected for exposure 
map effects and background. We adaptively bin the gas 
maps with an adaptive binning technique using weighted 
Voronoi tesselations (Diehl & Statlcr 2006; Cappellari & 
Copin 2003)3. 

For a subset of 36 galaxies with sufficient signal, we 
characterize the overall shape of the gas by deriving ellip- 
ticity and position angle profiles. A qualitative compari- 
son with optical DSS-2 i?-band images and a quantitative 
comparison with published optical surface photometry 
reveal little correlation between optical and X-ray prop- 
erties. In particular, we find no correlations between the 
optical and X-ray ellipticities at any radius. Even within 
one effective radius, where stars dominate the gravita- 
tional potential (e.g. Mamon & Lokas 2005; Humphrey 
et al. 2006), the two are completely uncorrelated. The 
second column of Table 1 lists mean gas ellipticities in 
the annulus between 0.6 and 0.9 J-band effective radii, 
from Paper I. These findings establish that the gas in el- 
liptical galaxies is, in general, not in precise hydrostatic 
equilibrium, at least at a level that precludes using the 
gas isophotes to constrain the shapes of dark matter halos 
(Buote & Canizares 1994; Buote et al. 2002). Instead of 
being hydrostatically calm, the gas morphology is almost 
always disturbed, even for rather X-ray faint galaxies. 

In this paper, we restrict our analysis to the subset of 
36 galaxies (Table 1) for which we derive ellipticity pro- 
files in Paper I. Following the technique used in Paper 
I to derive photon-flux calibrated images, we also pro- 
duce calibrated energy-flux images, which are necessary 
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Fig. 1. — Left panel: Gas X-ray surface brightness (in photons s"^ cm" ^ arcsec"^), eUipticity and position angle profile for NGC 4649. 
Solid linos indicate the Chebyshev polynomial fits. Right panel: Adaptively binned gas map of NGC 4649 with the fitted elliptical isophotes 
overlaid. 
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Fig. 2. — Gas surface brightness maps for three elliptical galaxies, showing the range in asymmetry index rj; Left; NGC 6482 {rj = 0.032); 
Center: NGC 4472 (77 = 0.077); Right: NGC 4374 (?7 = 0.478). Colors depict gas surface brightness levels with the color scale ranging from 
5 X 10-1° to 10-'5photonss^i cm ^ arcsec ^. The axes are labelled according to right ascension and declination (2000). 



to determine the central X-ray AGN Imirinosity (§2.3). 
We create mono-energetic exposm'e maps in steps of 7 in 
PI (~ 100 eV) and generate counts images in each indi- 
vidual PI channel (14.6 eV-wide). We then divide each 
counts image by the energetically closest exposure map 
to create a photon-flux-calibrated "slice," and multiply 
each slice by the appropriate energy of the PI channel. 
Finally, we sum all energy-weighted slices to produce a 
calibrated energy-flux image that automatically adjusts 
the effective exposure map to spatial changes in the spec- 
tral composition of the galaxies, due to radial temper- 
ature gradients and/or point source or AGN contribu- 
tions. This energy-flux calibrated image can be used to 
derive flux values, upper limits, and flux-calibrated sur- 
face brightness profiles even in situations with insufficient 
signal to perform spectral fits. 

2.2. Quantifying Asymmetry 
2.2.1. Definition of the Asymmetry Index rj 



Since many of the Chandra data sets lack the required 
signal and resolution to reliably identify individual asym- 
metric features in the gas surface brightness maps, we 
are forced to measure asymmetry statistically. We first 
construct a smooth, symmetric surface brightness model 
to subtract from the gas map. We base this model on 
the isophotal profiles computed in Paper I. These profiles 
give surface brightness (/x), eUipticity (ex)j and major 
axis position angle (PAx) as a function of mean radius 
r. We fit this surface brightness profile in logr — log/x 
space with a quadratic Chebyshev polynomial. For two 
galaxies with more complex surface brightness profiles 
(NGC 1316 and NGC 4636), we use a third order poly- 
nomial, and for one other object (NGC 4472), we resort 
to fourth order. To accommodate moderate changes in 
eUipticity and isophotal twists, we fit the eUipticity and 
position angle profiles in log r — ex and log r — PAx space 
with straight lines. From the fits, we compute a smooth, 
symmetric model for the gas distribution. The left panel 
of Figure 1 shows all three profile fits for NGC 4649. 
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TABLE 1 

Chandra X-ray gas morphology, AGN luminosity and temperature gradient 
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^Mean X-ray gas ellipticity ex between 0.6 — 0.9 Rj (J-band effective radius from 2MASS 
extended source catalog, Jarrett et al. 2000). 
''Asymmetry index r;, measured between Tmin and rmax (in kpc). 

•^X-ray AGN luminosity (in ergss~^) between 0.3 — 5keV, derived from and Sersic model 
plus point source fits to radial surface brightness profiles. 
''Mean logarithmic temperature gradient dlnT/d\nr evaluated at radii between 2 and 4 Rj 
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The adaptively binned gas map in the right panel has 
the elliptical isophotes of the adopted fits overlaid. Fi- 
nally, we bin the surface brightness model to match the 
binning structure of the adaptively binned gas map and 
subtract the binned smooth model to reveal small scale 
asymmetries in the gas distribution. 

To quantify the degree of asymmetry in this residual 
map, we define the asymmetry index rj in the following 
way: 




The asymmetry index is the sum over all pixels i of the 
squared relative deviations of the binned gas image Gi 
from the smooth model Mi, over and above the expected 
statistical deviations due to Poisson noise ctga^- Wc find 
that r] is an unbiased measure of asymmetry, independent 
of exposure time, radial fitting range, background level 
or signal-to-noise ratio. The errors on the rj values are 
estimated from a bootstrap analysis of 20 Monte-Carlo 
simulations. 

Figure 2 shows representative examples spanning the 
range of ry values present in our sample. NGC 6482 
(?7 < 0.032, left) shows no apparent signs of disturbance 
and is statistically consistent with a smooth, purely el- 
liptical model (Khosroshahi et al. 2004). In NGC 4472 
(77 = 0.077, center), one can detect low- level surface 
brightness depressions at several radii (Biller et al. 2004). 
NGC 4374 (right; Finogucnov & Jones 2001) has the 
highest r] value (rj = 0.48) in our sample. This sequence 
can serve as a guide for the meaning of 77. However, one 
should keep in mind that the examples in Figure 2 have 
been chosen to have high signal-to-noise ratios for display 
purposes. Most galaxies have significantly less signal and 
the asymmetric structures cannot be picked out by eye 
as easily. 

The third, fourth, and fifth columns of Table 1 list the 
values of rj and the range of radii from which they are 
derived. We have chosen to calculate rj from the full ra- 
dial range over which isophote fits are possible for each 
object, in order to take advantage of as much detected 
signal as possible. This choice unavoidably creates a cer- 
tain inhomogeneity in the sample since rj is not always 
being measured at the same physical scale. However, we 
have tested the sensitivity of all of the results of this pa- 
per to the size of the extraction annulus by separating 
into subsamples with small and large outer radii Tmax- In 
no case do we find any statistically significant difference 
between the subsamples. 

2.2.2. Simulations and Tests 

To test the behavior of the asymmetry index rj, we 
simulate images that qualitatively reproduce the variety 
of gas morphologies seen in the data. We start with a 
256 X 256 pixel realization of a Scrsic model with index 
n = 4 and a half-light radius of 50 pixels. We then add 
disturbances as relative surface brightness depressions 
surrounded by enhanced rim emission, mimicking the ap- 
pearance of AGN-induced "bubbles" , often observed in 

Equation (1) is equivalent to an area- weighted average over 
bins. We find empirically that an unweighted average gives too 
much emphasis to bright regions and small-scale structure. 



clusters of galaxies (e.g. McNamara et al. 2000). Fig- 
ure 3 shows a one-dimensional cut through our adopted 
circular bubble template profile, which is obtained by 
subtracting two e"^''/"'-' functions with slightly different 
a values, and adjusting their normalizations such that in- 
tegrating the 2-d profile yields a total deviation of zero. 
We define the "depression strength" as the maximum 
relative deviation at the bubble center, indicated by the 
— 100% mark in our template, while the dashed vertical 
lines in Figure 3 mark the "bubble size" . We then control 
the level of asymmetry in our simulations by varying the 
depression strength, bubble size, and the total number 
of bubbles. 

Wc simulate the disturbed gas images by randomly dis- 
tributing the bubbles over the entire Scrsic model image. 
Then, following our real data analysis, we construct a 
smooth model to subtract. For simplicity, we extract a 
radial surface brightness profile from circular annuli, and 
fit it with a Sersic model, instead of constructing a fully 
general elliptical model. We adaptively bin the image to 
a signal-to-noise per bin of 4, and apply the same bin- 
ning structure to the smooth model. After subtracting 
the binned model, we derive the asymmetry index from 
the residuals according to equation (1). We produce a 
suite of test simulations, keeping two parameters fixed 
while varying the third. For small numbers of simulated 
bubbles, the disturbed surface brightness models resem- 
ble galaxies in which one is able to reliably identify the 
location of individual bubbles. Although very large num- 
bers of simulated bubbles are unrealistic for real galaxies, 
we increase the number of bubbles beyond the physically 
motivated level to achieve more complex asymmetries in 
the morphology. 

Wc find that ?] is largely independent of the bubble size, 
but correlates positively with the number of simulated 
bubbles and their relative depression strength, as shown 
in the left and middle panel of Figure 4. We compute the 
expected trends by simulating 20 random spatial config- 
urations for each parameter set without adding Poisson 
noise, and computing the average of all rj values of these 
"perfect" data sets. We find that our asymmetry index 
values are fair representations of the expected values. 

The asymmetry index is constructed to account for the 
statistical scatter in the data and to be insensitive to dif- 
ferences in data quality. To verify this, we repeat our 
tests, keeping the bubble parameters fixed while varying 
the total galaxy luminosity, background level, cutoff ra- 
dius, and the exposure time of the observation. We find 
no correlations with any of these parameters. As an ex- 
ample, in the right panel of Figure 4 we show one test in 
which we vary the exposure time. One can see that 77 is 
not affected by the increasing resolution that is a result 
of the rising total number of counts. We conclude that 
77 is an unbiased measure of the asymmetry present in 
the gas maps, and unaffected by non-morphological data 
properties. 

2.3. X-ray AGN Luminosities 

We identify the location of the central AGN in the X- 
ray image by overlaying the high-accuracy center from 
the 2MASS extended source catalogue (Jarrett et al. 
2000). The 2MASS position is derived from the lumi- 
nosity weighted center of the co-added K, J and H band 
images and has a typical Icr-unccrtainty of only 0.3". We 
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Fig. 3. — Onc-dimonsional cut through a simulated "bubble de- 
pression" , demonstrating the fractional depression strength of the 
bubble with respect to its radial extent. The bubble has negative 
surface brightness deviations within the bubble radius RBubble> a^nd 
slightly enhanced rims around it. The profile shape is chosen such 
that the integrated deviations integrate to over area. 
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Fig. 4. — Tests of the asymmetry index using models with sim- 
ulated bubble-like depressions (see Fig. 3). The asymmetry index 
77 is sensitive to the number of simulated bubbles (left), and the 
relative strength of the depression (middle), but not to the effec- 
tive exposure time, characterized by the total number of counts 
(right). The dashed lines mark the expected behavior for perfect 
(noise-free) data. 

combine this uncertainty with Chandra^s pointing accu- 
racy^, and manually assign the AGN to the brightest X- 
ray point source in the 0.3 — 5.0 keV counts image within 
a ~ 3 pixel radius of the 2MASS center. If multiple 
source candidates are within this field, we take the point 
source closest to the centroid of the diffuse emission. In 
cases where soft gas emission completely dominates the 
central region, we revert to the 2 — 5 keV band for iden- 
tification instead. This technique allows a rather reliable 
AGN identification, with only few ambiguous cases. 

Unfortunately, the location of the AGN coincides with 
the peaks of both the diffuse gas emission and the stel- 
lar point source component. While this sometimes helps 
in pinpointing the location of the AGN, it makes it diffi- 
cult to determine the AGN flux. The standard technique 
is to extract a spectrum from the central region and to 
attempt to spectrally disentangle the AGN component 
from the contaminating diffuse gas and low-mass X-ray 
binaries (LMXBs) (e.g. Kim & Fabbiano 2003). How- 
ever, the weak AGN in normal elliptical galaxies often 
yield insufhcicnt counts to constrain all parameters in 
this complex spectral fit. This fit is additionally compli- 
cated by the fact that the AGN and unresolved LMXB 
components are both well-described by featureless power- 
laws, and can sometimes become virtually indistinguish- 
able. For many observations in our sample, the AGN 
contributes no more than a few counts, making a reli- 
able spectral fit practically impossible. 

^ http: / / cxc.harvard.edu / mta/ASPECT/abs_point.html 



Instead, we decide to spatially disentangle the X-ray 
AGN emission from its complex background using the 
energy-flux calibrated images of diffuse emission. The 
energy weighting makes it easier to distinguish the hard 
AGN component. We bin the energy-flux image into cir- 
cular annuli centered on the AGN position and adap- 
tively change the annular bin sizes to enforce a minimum 
signal-to-noise requirement of 2. As this would result in 
very large numbers of bins at larger radii, we also require 
the annular width to be larger than 10% of the mean bin 
radius to ensure proper azimuthal averaging within a bin, 
generally resulting in a logarithmic binning at large radii. 

We fit the inner 32 arcsec of the radial surface bright- 
ness profile with a two-component model, one represent- 
ing the point source and one the diffuse emission. The 
point source profile is obtained from a normalized mono- 
energetic point-spread function (PSF) centered on the 
position of the central AGN and computed at its mean 
photon energy. We bin the PSF image to the same cir- 
cular binning, to give a radial point source model that 
matches our binning structure with only the normaliza- 
tion flux as a free parameter. To represent the diffuse 
emission we use a Scrsic or a /? model, on top of a uni- 
form background. Figure 5 shows the inner 10 arcsec of 
two galaxies with a statistically signiflcant AGN detec- 
tion (NGC 4261, left) and a non-detection (NGC 6482, 
right). The dashed lines indicate the diffuse component 
from the Sersic (grey) and f3 model fit (black), while the 
solid lines show the same fits with the additional AGN 
component. All fits are inspected by eye and obviously 
unconverged fits are repeated over a different radial ex- 
tent until a satisfactory fit is obtained. 

Even though both models generally produce statisti- 
cally equally good fits, the (3 model fits systematically 
yield higher fluxes for the AGN, as the /?-proflle is flatter 
at small radii. We adopt the average of the AGN fluxes 
yielded by the (3 and Scrsic model fits as the best value. 
The difference between the two individual fits and the 
average value is assumed to represent a 3a systematic 
error and is combined with the statistical errors yielded 
by each profile fitting procedure. Finally, we correct the 
AGN fiuxes for Galactic absorption with the same correc- 
tion factor that is used to correct the total gas flux, com- 
puted from integrating the best spectral fit for the gas 
emission with and without the Galactic absorption fac- 
tor. This correction may be a slight overestimate, since 
the AGN emission is generally harder and less strongly 
affected by absorption; on the other hand the column 
density is also likely to be higher due to intrinsic absorp- 
tion inside the galaxy. This correction is generally rather 
small, well below our statistical and systematical errors, 
and does not affect the computed X-ray AGN luminosity 
significantly. 

The final AGN luminosities Lx.agn are reported in 
the sixth column of Table 1, together with the combined 
statistical and systematic errors. Flux values that are de- 
tected below the la confidence limit are reported as up- 
per limits. For objects with insufficient signal to produce 
a surface brightness profile, we sum up the energy-flux 
calibrated image within the region containing 99% of the 
PSF fiux and correct for a fiat background extracted from 
an adjacent annulus extending out to 10 arcsec. Since we 
have no way to distinguish between a central peak in gas 
emission, LMXBs, or an actual AGN for these galaxies. 
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we quote the 3cr upper bound of the summed flux as an 
upper limit on their AGN luminosity. 

2.4. Radio Luminosities 

We use the NRAO VLA Sky Survey (NVSS; Condon 
et al. 1998) to derive homogeneous 20 cm radio contin- 
uum luminosities for our sample. As we are interested in 
the AGN's impact on the ISM, we sum up all associated 
radio sources within 3 J-band effective radii from the 
2MASS catalog (Jarrett et al. 2000) radu to derive the 
NVSS radio luminosity iNVSSi thus including contribu- 
tions from extended structures. This results in general 
in an extraction radius of approximately 1.5' for most 
of the sample (see Paper I for a complete list of optical 
radii). The formal average NVSS source density for the 
full NVSS catalog is about 50 sources per square degree 
(Condon et al. 1998), resulting in a probability to include 
unrelated background sources within 1.5' of only 9%. For 
63% of our sample, this probability lies below 10%, for 
91% it lies below 20%. To further reduce the chance of 
including background sources, we inspect all NVSS im- 
ages by eye and compare radio features with X-ray and 
optical images to find possible counterparts. We manu- 
ally remove obviously non-associated background sources 
from the source list. 

The resulting luminosities are listed in the second col- 
umn of Table 2. For the majority of our sample, the 
limited spatial resolution of the NVSS survey precludes 
reliably distinguishing between extended and point-like 
radio sources. Only where unambiguously possible, we 
determine the NVSS position angle by manually measur- 
ing the orientation of the major axis. 

Twenty-three galaxies in our sample have been ob- 
served in the Faint Images of the Radio Sky at Twenty- 
Centimeters (FIRST) VLA survey (Becker et al. 1995), 
which has a much higher spatial resolution than NVSS, 
with a similar detection limit. We extract FIRST ra- 
dio luminosities from the inner 3 optical effective radii 
to match the NVSS extraction region. A comparison be- 
tween the two radio surveys indicates that both are yield- 
ing consistent values. However, the higher spatial reso- 
lution of the FIRST survey gives us the opportunity to 
extract fluxes from the central point source alone, largely 
excluding contributions from the radio jets. Thus, in- 
stead of using a large extraction radius, we sum the flux 
within a smaller 30 arcsec radius to determine the FIRST 
luminosity ipiRST, listed in the third column of Table 2. 
These values are systematically smaller than the NVSS 
fluxes due to the smaller extraction region, but still very 
well correlated. This indicates that our larger NVSS sam- 
ple is not heavily contaminated by background objects, 
but that the extra flux is most likely associated with ex- 
tended radio structure intrinsic to the galaxy. 

2.5. Correlation between Radio and X-ray AGN 
Luminosities 

Figure 6 compares the AGN X-ray luminosity with 
both the NVSS (left) and FIRST (right) 20 cm contin- 
uum radio powers. Both plots show a clear correlation 
between radio and X-ray luminosities, spanning almost 
six orders of magnitude. The slope is consistent with a 
purely linear relation between both luminosities, as indi- 
cated by the solid lines. We analyze all of the correlations 



TABLE 2 

Radio luminosities and galaxy environment 



Name -T-NVBs"^ ipiRST'" logio P2MASs'' 



''20 cm continuum radio luminosity (in ergss ^ Hz ^) from NVSS 
within 3Rj, and FIRST within 30". 

''Projected local galaxy density (in Mpc"^)^ derived from the 
2MASS extended source catalog (Jarrett et al. 2000). Reported er- 
rors only include statistical errors, not systematic errors. 
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Fig. 5. — Radial X-ray surface brightness profiles for the inner 6" of NGC 4261 (left) and NGC 6482 (right). Thick solid lines indicate 
Sersic model fits, thin solid lines show /3 model fits, including a PSF model for the central AGN. Dashed lines show the same model fits 
without the AGN component. NGC 4261's AGN is detected at a Tcr level, NGC 6482 profile reveals no AGN signature. For NGC 6482, 
the AGN flux yielded by the Sersic model fit is so small that the thick dashed line is invisible, as it is covered by the thick solid line. 



in this paper using an algorithm ("bandfit"; see Ap- 
pendix A) that models the distribution of data points as a 
linear band with finite Gaussian intrinsic width. Bandfit 
can handle data with errors in either or both variables, 
as well as censored data (upper or lower limits). The 
bandfit analysis for the iNVSS^^x,AGN relation puts 
the chance for the null hypothesis of no correlation to 
less than 10~^^% Sa). This correlation is not sim- 
ply a consequence of plotting distance vs. distance; it is 
equally significant if the distance dependencies of both 
parameters are removed. The best fit correlation can be 
described as 

logLx,AGN = 0.62(±0.10) logLNvss + 19.11(±0.12). 

Similar correlations have been found in the past for 
AGN (e.g. Hardcastle & Worrall 1999; Brinkmann et al. 
2000, and references therein), with slopes ranging from 
~ 0.5 for radio-loud quasars to > 1 for radio-quiet 
quasars and radio galaxies. One must remember, how- 
ever, that these results were obtained using total X-ray 
luminosities, not nuclear lumninosities as in our case. 
Also, previous samples have consisted mostly of objects 
in which the AGN dominate the X-ray emission, in con- 
trast to our sample where the AGN is outweighed by the 
gas and LMXB contributions. If we replace Lx,agn with 
the total gas luminosity ix,gasj we find a steeper slope, 
similar to previous results for radio galaxies. 

Hardcastle & Worrall (1999) have argued that the ex- 
istence of an L^r^^radio Correlation indicates that both 
X-ray and radio emission originate in a strongly Doppler 
boosted jet. If this were the case, then measured X-ray 
or radio fiuxes might refiect primarily the orientation of 
the jet relative to the line of sight, rather than the intrin- 
sic power of the AGN. In this picture, the shallower slope 
of the correlation for quasars could be explainable by a 
combination of beamed and isotropic X-ray emission, im- 
plying that X-rays would, at least in part, measure AGN 
power. Whether the shallow slope we find in the -Lnvss^ 
^x.AGN relation indicates a similarity with the quasars is 
unclear. It is possible that our slope could be affected at 
the low-luminosity end by occasional misidentification of 
individual LMXBs in the galaxy centers as AGN. How- 



ever, if Lnvss or Lx.AGN were only reflecting the jet 
orientation, we would not expect either to be correlated 
with gas morphology, as we do see. 

A survey of radio emission from normal elliptical galax- 
ies by Wrobel & Heeschcn (1991) suggests that the ma- 
jority of their radio fiux is due to AGN activity rather 
than star formation. The fact that nuclear X-ray and ra- 
dio luminosities show a tight relationship supports this 
conclusion for our sample, as well as indicating that we 
have properly identified and isolated the central AGN in 
the X-ray images. We conclude that both the radio and 
the X-ray AGN luminosities arc reliable tracers of the 
central AGN activity in our sample, most likely due to 
jets powered by gas accreting onto a supermassive black 
hole. 

2.6. The Galaxy Environment 

To quantify the likelihood of interactions with neigh- 
boring galaxies, we measure the projected galaxy den- 
sity in the vicinity of our sample galaxies by counting 
the number of neighbors listed in the 2MASS extended 
source catalog (Jarrett et al. 2000). We identify all galax- 
ies within a projected radius of 100 kpc at the distance 
of the target as "neighbors." We restrict this distance 
to a maximum angular size of 15 arcmin, to minimize 
contamination by random foreground or background ob- 
jects. Despite the smaller physical extraction radius for 
nearby objects, the number of neighbors is still sufficient 
to derive a galaxy density, as the incompleteness limit 
drops significantly, and lower-luminosity neighbors can 
be identified. Because the 2MASS catalog is flux-limited, 
we have to cope with a distance-dependent completeness 
limit. We use the K-hand Schechter luminosity func- 
tion from Kochanck et al. (2001) to correct for incom- 
pleteness. We integrate the luminosity function down 
to the completeness limit at the object's distance and 
divide by the integral of the luminosity function down 
to a reference luminosity, which we arbitrarily chose as 
Mk = —11.5. This way we can calculate the fraction of 
all galaxies above the reference luminosity that 2MASS is 
able to detect. We find that, after correcting for incom- 
pleteness, our galaxy densities P2MASS a-rc independent 
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Fig. 6. — X-ray AGN luminosity vs. 20 cm continuum radio luminosity with data taken from NVSS within 3 optical radii (left) and from 
FIRST (right) within a 30" radius. Arrows indicate 3cr upper limits. The slope of the relationship between radio and X-ray luminosity is 
close to linear, as indicated by the solid lines. The dashed line in the left panel shows the bandf it best fit relation, which has a slope of 
0.62 ± 0.10. Thus, X-ray and radio luminosities are each measures of AGN activity. 



of distance. We list the values for p2MASs in Table 2 
in units of galaxies per Mpc. The errors include only 
statistical errors due to the galaxy counting statistics; 
deriving systematic errors due to background objects or 
uncertainties in the galaxy luminosity function are be- 
yond the scope of this paper. 

Ideally, one would also like to restrict the list of neigh- 
bors in velocity space, to remove galaxies that are simply 
aligned along the line of sight. Such a three-dimensional 
restriction requires an additional unbiased source of ra- 
dial velocities. However, available radial velocities are 
heavily biased toward "more interesting" regions of the 
sky and better studied objects. Thus, we refrain from 
using these velocities and use the projected galaxy den- 
sity P2MASS instead, a measure similar to the TuUy den- 
sity parameter (TuUy 1988), but well-defined for our en- 
tire sample. All conclusions based on P2MASS £^re repro- 
ducible with the TuUy parameter instead. 

The hot ISM morphology may also be affected by hy- 
drodynamic interactions with a surrounding intra-cluster 
medium (ICM) or intra-group medium (IGM). One ex- 
pects a transition to an IGM or ICM at large radii to be 
accompanied by an outward rise in gas temperature, re- 
flecting the higher virial temperature of the group or clus- 
ter. To assess the presence of a hotter ambient medium, 
we compute radial temperature profiles for our sample. 
The X-ray counts image of each galaxy is divided into el- 
liptical annuli, according to the X-ray ellipticity profiles 
computed in Paper I, except for those with insufhcient 
signal to fit ellipses, for which we revert to circular an- 
nuli. We then extract a source and background spectrum 
for each annulus and fit them with a two-component 
model in the CIAO analysis package Sherpa. The first 
component consists of an APEC ^ plasma model to rep- 
resent the hot gas emission. A quantitative compari- 



son with its better known predecessor, the Mekal model, 
shows nearly identical results. We fix the gas metal- 
licity at the solar abundance value. Unresolved point 
sources are represented by a power-law model with the 
power law index fixed at 1.6. This "universal" spectral 
model is an adequate representation of the emission of 
low-luminosity low-mass X-ray binaries, as demonstrated 
in Paper I and determined independently by Irwin et al. 
(2003). We also add a multiplicative absorption compo- 
nent, for which we fix the hydrogen column density to 
the Galactic value, evaluated at the target position with 
the CIAO tool ColderJ . We repeat our spectral analy- 
sis for a few objects with the gas abundance as a free 
parameter, and find that our choice to fix them to the 
solar value does not affect the fitted temperature. Since 
the mctallicity is poorly constrained by the fits in low 
signal-to-noise systems, we fix the metallicity for all of 
our galaxies, in order not to introduce systematic dif- 
ferences in the analysis. We then fit the temperature 
profile by a power law, and derive a mean logarithmic 
temperature gradient between 2 and 4 J-band effective 
radii (i?j): a24 = dhiT/dlnR\2-4Rj. Values of Q!24 are 
listed in the last column of Table 1. 

One can consider the galaxy density /32MASS ^-s an in- 
dicator of the likelihood of gravitational galaxy interac- 
tions through close encounters or mergers, similar to the 
TuUy galaxy density parameter. The outer temperature 
gradient, on the other hand, is a measure of the presence 
of hot ambient gas and thus a proxy for the likelihood of 
hydrodynamic interactions with this gas. A much more 
thorough analysis of the temperature profiles present in 
our sample will be the main focus of Paper III. 

3. RESULTS 
3.1. X-ray Ellipticity and Asymmetry 
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Fig. 7. — The X-ray asymmetry index as a function of X- 
ray ellipticity ex- The observed positive trend suggests a common 
underlying cause for both the large measured X-ray cUipticities and 
the asymmetries. 

Figure 7 shows a comparison between the mean X- 
ray eUipticities between 0.6 — 0.9 optical radii and the 
asymmetry index rj. There is a trend for flatter galaxies 
to be also more asymmetric, with a bandf it correlation 
strength of 99.72%. There is a clear deficit of galax- 
ies with flattened X-ray emission, but small asymmetry, 
which suggests a common underlying cause for both. The 
ellipticity may best be interpreted as simply a mean mea- 
sure of the m = 2 Fourier amplitude of the departure 
from circular symmetry. This is consistent with our re- 
sults from Paper I, where we find that the gas eUiptici- 
ties are uncorrelated with the eUipticities of the starlight, 
and thus are not simply determined by the shape of the 
smooth underlying gravitational potential. 

3.1.1. AGN Influence on Gas Morphology 

Figure 8 shows the dependence of the asymmetry index 
77 on the two independent measures of AGN activity: the 
X-ray luminosity ix,AGN of the central point source (left 
panel) and the NVSS radio power Lnvss > integrated over 
three optical radii (right panel). Both AGN properties 
are correlated with rj, in the sense that galaxies hosting 
stronger AGN have stronger morphological asymmetries 
in their hot ISM. Analysis with bandf it rejects the null 
hypothesis of no correlation between Lnvss a-nd ry at the 
98.7% confidence level (sec Table 3). The same analysis 
for Lx,AGN and rj yields a slightly lower confidence value 
of 97.1%. If we regard Lx,agn and Lnyss as indepen- 
dent measures of the same phenomenon (§2.5), we can 
compute the joint probability of the null hypothesis of 
no correlation between general AGN activity and asym- 
metry, which is then rejected at the 99.96% level. Sim- 
ilar results are obtained from bootstrapping and from a 
Spearman rank analysis ignoring the upper limits. This 
AGN-asymmetry connection is strong evidence that the 
central AGN is generally responsible for disturbing the 
hot gas morphology even in normal elliptical galaxies and 
provides further support for the importance of AGN feed- 
back. 



TABLE 3 

Significances for correlations with tj. 



Second Variable PnuU 



logLNVSS & log Lx, AGN (jointly) 0.04% 

024 0.2% 

logLNVSS 1-3% 

logLx.AGN 2.9% 

logP2MASS 24.9% 

logPTuiiy 45.1% 



A comparison with our simulations in §2.2.2 suggests 
two possible interpretations for the increase of asymme- 
try with AGN luminosity: the number of inflated bub- 
bles, or the strength of the created depressions. The 
asymmetry index is equally sensitive to both and cannot 
distinguish between them. However, a qualitative inves- 
tigation of the gas morphologies (see Paper I) suggests 
that depression strengths are more important. The ex- 
amples in Figure 2 already show that galaxies with larger 
ry values tend to have stronger surface brightness depres- 
sions, possibly caused by the stronger interaction of more 
powerful radio sources with the ISM. 

3.2. Environmental Influence on Gas Morphology 

The probability of interactions with neighboring galax- 
ies through mergers, tidal interactions or close encoun- 
ters naturally increases with the local number density of 
galaxies. Thus, if interactions are important in influenc- 
ing the gas morphology of elliptical galaxies, we would 
expect a correlation between morphological parameters 
and the projected galaxy density /32MASS- However, Fig- 
ure 9 shows that neither the gas ellipticity (upper panel) , 
nor the asymmetry index (lower panel) shows any trend 
with this measure of environment. A bandf it analysis 
yields no statistically significant correlation in either case 
(Table 3). Bootstrapping and a Spearman rank analysis 
give consistent results. 

To ensure that our conclusions are independent of our 
definition of galaxy density, we repeat the analysis with 
other environmental measures of galaxy density: the 
TuUy density parameter pTuiiy (TuUy 1988), the distance 
to the 10th closest neighbor in the 2MASS extended 
source catalog (Jarrett et al. 2000), and the number of 
galaxies associated with a galaxy group, as listed in the 
Lyon Group of Galaxies (LGG; Garcia 1993) catalog. All 
of these measures of environment yield the same result: 
the galaxy density is not correlated with the observed 
asymmetries or large eUipticities of the hot gas. Thus, 
the X-ray gas morphology is not mainly driven by in- 
teractions with companions, although they may still be 
important for some individual objects. 

In the absence of galaxy-galaxy interactions, hydro- 
dynamic interaction with a hot ambient medium is still 
possible, through ram pressure, shocks, pressure confine- 
ment or hydrodynamic instabilities at the contact inter- 
face. Two galaxies in our sample are already known to 
be affected by ram-pressure: NGC 4472 (Biller et al. 
2004) and NGC 1404 (Machacek et al. 2004). However, 
in both cases, the radius within which we compute the 
asymmetry index is smaller than the radius at which the 
interaction with the ambient medium becomes obvious. 
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Fig. 8. — Asymmetry index r) as a. function of X-ray AGN luminosity Lx.AGN (left) and 20 cm radio continuum power LnvSS taken from 
NVSS (right). Both plots indicate that the central AGN luminosity is positively correlated with the observed asymmetry in the gas, with 
the correlation extending all the way down to the weakest AGN. 
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Fig. 9. — X-ray ellipticity ex (top) and asymmetry index ri 
(bottom) as a function of projected galaxy density P2MASS- Both 
plots indicate that environment is not the driving factor causing 
the ellipticities and asymmetries. 

To assess the sample in general, we use the outer tem- 
perature gradient 0:24 as a proxy for a hot ICM, IGM, or 
circumgalactic gas. A larger gradient suggests a larger 
ambient gas pressure. Figure 10 shows that asymmetry 
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Fig. 10. — Asymmetry index 77 as a function of the outer temper- 
ature gradient a24 1 an indicator for the presence of hotter ambient 
gas. The correlation is such that galaxies in a hot gas environment 
arc more asymmetric. Note that the sample includes both rising 
(a24 > 0, hotter ambient gas present) and falling (024 < 0, no 
ambient gas present) temperature profiles. 

indeed increases as a function of Q!24. A bandf it analy- 
sis rejects the null hypothesis at the 99.8% level (Table 
3), with similar results obtained from bootstrapping and 
Spearman rank correlations. 

4. DISCUSSION 

4.1. Causes of Disturbances in the Gas Morphology 

Our comparison of optical and X-ray morphologies in 
Paper I revealed that even well within one effective ra- 
dius, where the gravitational potential is dominated by 
the stellar component (Mamon & Lokas 2005; Humphrey 
et al. 2006), there is no correlation between stellar and 
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gas ellipticities, very much at odds with predictions from 
hydrostatic equihbrium. In fact, many X-ray eUipticities 
even exceed those of the starhght, quite the opposite of 
what is expected from purely gravitational effects. Wc 
also noted in Paper I the prevalence of asymmetries in 
the hot gas, which we have now statistically quantified in 
this paper with the asymmetry index rj. We find that the 
large observed ellipticities and the amount of asymme- 
try in the gas are correlated, pointing toward a common 
cause that dictates the overall X-ray morphology. 

In this Paper, we have shown that feedback from the 
central AGN (Figure 8) and hydrodynamic interactions 
with the ambient medium (Figure 10) are the dominant 
causes of these asymmetries. Our analysis of the AGN- 
morphology correlation shows that this phenomenon is 
not restricted to powerful radio sources, but forms a con- 
tinuous sequence down to the weakest AGN. The corre- 
lation with the outer temperature gradient a24, on the 
other hand, suggests that the presence of a hot am- 
bient medium has a similarly strong effect on the gas 
morphology. Figure 11 shows that both correlations to- 
gether dictate the disturbance in the gas. In Figure 11a 
we have colored the residuals in the a24 — 77 correla- 
tion according to AGN luminosity, with red indicating 
low (Xnvss < lO'^^crgs^^ Hz"^), and green high, ra- 
dio power. In Figure lib we color the residuals in the 
-^^Nvss— correlation according to a24, with red now indi- 
cating a negative and green a positive temperature gradi- 
ent. (Black points indicate objects for which these values 
are unavailable.) In each case, the residuals are clearly 
correlated with the third quantity. Replacing Lnvss with 
the X-ray AGN luminosities £x,agn yields comparable 
results. Note that the scatter about the mean correla- 
tions is about the same in each case, indicating that AGN 
and ambient medium have comparable effects on the gas 
morphology. 

4.2. Implications for AGN Feedback in Elliptical 
Galaxies 

The great majority of galaxies in our sample are not 
conspicuous AGN; most harbor only weak, unresolved ra- 
dio sources. Nevertheless, we find a correlation between 
hot gas asymmetries and the radio and X-ray AGN lu- 
minosities, down to the lowest luminosities detectable by 
NVSS, FIRST, or Chandra. These results solidify the 
AGN as an influential feedback mechanism for normal 
elliptical galaxies in general, even where they do not ob- 
viously harbor powerful nuclear sources. 

Intermittent AGN activity may stir the hot gas, push- 
ing it far enough out of equilibrium that information 
about the shape of the underlying potential is lost. That 
we observe morphological disturbances in nearly all ellip- 
tical galaxies, and that some objects even reveal signa- 
tures of multiple outbursts (e.g. O'SuUivan et al. 2005), 
suggest that the duty cycle of the AGN is shorter than 
or comparable to the sound crossing time. Thus, the 
hot gas is continually disturbed, without having time to 
resettle into equilibrium. 

Despite the lack of correlation between optical and X- 
ray ellipticities (Paper I), we do note a weak tendency 
for the hot gas isophotal major axes to be aligned with 
those of the starlight. To quantify this trend, we extract 
a subsample of 15 galaxies that have X-ray ellipticities 
ex > 0.2 in order to have well-defined position angles 



(PA). The PA differences between the gas and the stel- 
lar component, in 15°-wide bins, are shown as the grey- 
shaded histogram in Figure 12. A Kolmogorov-Smirnov 
(KS) test puts the probability for a chance correlation at 
only 0.05%. 

While this alignment may refiect the orientation of the 
underlying gravitational potential, it may also be the re- 
sult of interactions with the central radio source. Our 
sample exhibits an anti-correlation between radio and op- 
tical major axes, which is shown as the solid histogram 
in Figure 12. This controversial phenomenon has been 
observed early on by Palimaka et al. (1979), but has 
not been confirmed in subsequent studies (Sansom et al. 
1987). Nevertheless, a KS-test yields only a 1.1% proba- 
bility for the null hypothesis of uncorrclated orientations 
in our small sample of 11 galaxies with both available 
radio and optical major axes. 

In case of a relatively stable radio jet axis, one would 
expect the jet to preferentially evacuate gas along the 
axis, creating a misalignment between radio and X- 
ray gas emission. We do observe a very slight anti- 
correlation, shown in the dashed histogram in Figure 12, 
but only 12 galaxies have both reliable X-ray and radio 
position angles. The large error bars on both radio and 
X-ray position angles make it difficult to ascertain the 
significance of the correlation, which a KS-test puts nom- 
inally at 90%. Extracting the X-ray position angles at a 
larger radius increases the significance; but the detection 
of this misalignment is marginal and awaits confirmation 
with a larger sample and deeper radio data. With our 
current small sample, a two-sided KS test yields a 22.3% 
probability that the radio-X-ray anticorrelation and the 
X-ray-optical alignment are drawn from the same distri- 
bution, which rises to 80.6% at a larger extraction radius. 
Thus, the data are consistent with the alignment being 
caused by the effects of the radio source on the hot ISM. 
The competing effect of the interaction with the ambi- 
ent medium may also work against a proper detection 
of this misalignment. However, the position angle may 
also represent the only information that is still inferable 
about the underlying gravitational potential, namely its 
orientation. 

4.3. Implications for Interactions with the Intergalactic 

Medium 

We find that the observed ellipticities and asymmetries 
in the hot gas are not mainly a consequence of inter- 
actions with neighboring galaxies. On the other hand, 
interactions with the ambient medium seem to be as im- 
portant as the interplay with the central radio source. 
We expect the impact of an external ambient medium 
to increase with distance from the center. However, our 
analysis is confined to a region within Chandra^s limited 
field of view, restricted further by our need for sufficient 
signal to fit isophotes. Thus, our asymmetry index is in- 
sensitive to structure outside a radius which depends on 
the average surface brightness of the object and the ex- 
posure time of the observation. To verify that the — 
correlation is not simply a consequence of a variable ex- 
traction radius, we have checked for a correlation be- 
tween this radius and 77 and we find none. 

The majority of X-ray halos associated with elliptical 
galaxies extend well beyond the outer radii in Tabic 1, 
as was already demonstrated with observations with the 
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Fig. 11. — Similar to Figures 10 and 8. Left: correlation between rj and the outer temperature gradient 024- Red colors indicate 
low-luminosity (Lnvss < 10''^ ergs~^ Hz^^), and green colors higher luminosity (Lnvss > lO'''^ erg s~^ Hz~^) AGN. Right: correlation 
between r] and radio luminosity LnvSS- Red points indicate a negative, and green a positive, temperature gradient 024. In each case the 
residuals from the mean relation are correlated with the third variable. On can consider the more fundamental correlation to be a 2-d plane 
in log LNVSS-cr24-log r? space. 

tion of AGN in richer environment are switched "on" 
compared to isolated galaxies. In this scenario, the cor- 
relation between the outer temperature gradient and the 
gas asymmetry would be a result of an increased duty 
cycle for AGN outbursts, rather than hydrodynamic in- 
teractions with the ambient medium. In addition, the 
circumgalactic gas reservoir also acts as a pressure con- 
finement, which would prevent radio jets and inflated 
cavities from simply leaving the system and provides a 
background contrast to actually identify asymmetries at 
larger radii. 
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Fig. 12. — Histogram of absolute differences between optical 
(PAopt), radio (NVSS 20cm continuum, PAimvss)i ^i^d X-ray gas 
(PAx) position angles. The grey shaded histogram with dotted 
contours shows a tendency for galaxies to have X-ray and optical 
isophotes aligned. The solid line shows an anti-correlation between 
optical and radio position angles, and the dashed line a slight anti- 
correlation between radio and X-ray position angles, suggesting 
that the radio source is responsible for the X-ray-optical alignment. 



Einstein and iZ 05^ T satellites (e.g. Forman et al. 1985; 
O 'Sullivan et al. 2003). Interactions with the ambient 
medium may very well have an even stronger influence 
on the gas structure at larger radii. At this point we 
are not yet able to distinguish the effects of AGN from 
those of environment based on morphology alone. Nev- 
ertheless it is curious to note that we do not observe a 
correlation between the X-ray ellipticity and the outer 
temperature gradient, as one might expect. Perhaps a 
different deflnition of statistical asymmetry that is sen- 
sitive to "lopsidedness" may be useful in this regard. 

Alternatively, the existence of a positive outer temper- 
ature gradient could also be interpreted as the presence 
of an external fuel reservoir for the central AGN. Best 
(2004) finds a correlation between radio-loud AGN ac- 
tivity and environment, in the sense that a larger frac- 



5. CONCLUSIONS 

We have extended our earlier work from Paper I on nor- 
mal elliptical galaxies, which showed that optical and X- 
ray gas ellipticities are uncorrelated in the regions where 
stellar mass dominates, contrary to what is expected 
from gas in hydrostatic equilibrium. Instead, the hot gas 
is generally disturbed. To elucidate the nature of these 
asymmetries, we have introduced the asymmetry index rj, 
which measures the statistical deviation from a symmet- 
ric surface brightness model. This allows a systematic 
study of asymmetry even in relatively low-luminosity ob- 
jects, in which asymmetric features are barely resolved. 

We find a strong correlation between the asymmetry 
index and two independent measures of AGN activity: 
the radio continuum power at 20 cm from NVSS (Con- 
don et al. 1998) and the X-ray AGN luminosity extracted 
from Chandra data. The observed AGN-asymmetry cor- 
relation persists all the way down to the weakest AGN, 
where the NVSS survey reaches its detection limit. This 
is quite surprising, since these objects generally lack ex- 
tended jet signatures in their radio images and are mostly 
detected as weak central point sources, if at all. We 
also find the asymmetry index to be correlated with the 
outer temperature gradient 024, a proxy for the pres- 
ence of a hot intragroup or intracluster medium. The 
strength of this correlation is comparable to that of the 
AGN-asymmetry correlation, indicating that hydrody- 
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namic interaction with the ambient medium is compa- 
rably important to AGN in determining the X-ray mor- 
phology of normal ellipticals. However, we find no such 
correlation with the density of neighboring galaxies, and 
no evidence that galaxy-galaxy interactions play a sig- 
nificant role in shaping the hot ISM in these systems. 
Alternatively, the presence of an outer temperature gra- 
dient may also indicate the availability of a fuel reservoir 
for the central AGN. In this case, the increased amount 
of asymmetries in the X-ray gas could indicate a greater 
duty cycle for AGN outbursts in these objects. 

The emerging picture is consistent with the AGN per- 
sistently stirring up the interstellar medium through in- 
termittent outbursts, and strengthens the case for the 
AGN to be at least partly responsible for offsetting cool- 
ing in normal elliptical galaxies. We will address the im- 
pact of the central AGN on temperature and entropy pro- 
files in detail in Paper 111 of this series (Diehl & Statler 



2008). 
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APPENDIX 

CORRELATION ANALYSIS WITH BANDFIT 

The bandfit algorithm fits a model to two-dimensional {x,y) data that are, or are suspected of being, linearly 
correlated. The model is a linear band, described by an orientation angle 4> relative to the x axis, a perpendicular 
offset from the origin zq, and an intrinsic Gaussian perpendicular width a-mt- The likelihood of any observed datum 
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{xi,yi) is a function of the perpendicular distance to the midline of the band, and is straightforwardly calculated 
by convolving the band with the error ellipse for the datum. The algorithm is thus similar to line fitting using 
perpendicular residual minimization, which is not exactly the same as classical linear regression. 

Bandfit is written to handle data with Gaussian errors in either or both variables, including correlated errors. 
Censored data (upper or lower limits in either or both variables) are also easily handled, and contribute to the 
likelihood via an integral over the convolved band up to their location (i.e., the likelihood is proportional to the erf 
function evaluated at the censored point). 

The model parameters are found by maximizing the likelihood using the well known simplex ("amoeba") method. 
Errors in the fitted parameters are returned in the form of a covariance matrix. 

To assess a correlation, the data are first shifted to zero mean to minimize the covariance between <j) and zq. The 
square root of the (jycj) element of the covariance matrix can then be interpreted as the standard error in </>. Tests with 
simulated data sets show that integrating over this Gaussian provides accurate confidence limits on the non-zero slope 
of the correlation. 

Bandfit is implemented in IDL, and will be made available through the IDL astro library after a thorough workout 
by volunteer beta-testers. 



